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The 3d electron spinel compound L1V2O4 exhibits heavy lermion behaviour below 30K which 
is related to antiferromagnetic spin fluctuations strongly enhanced in an extended region of 
momentum space. This mechanism explains enhanced thermodynamic quantities and nearly critical 
NMR relaxation in the framework of the selfconsistent renormalization (SCR) theory. Here we 
show that the low-T Fermi liquid behaviour of the resistivity and a deviation from this behavior for 
higher T may also be understood within that context. We calculate the temperature dependence of 
the electrical resistivity p(T) assuming that two basic mechanisms of the quasiparticle scattering, 
resulting from impurities and spin-fluctuations, operate simultaneously at low temperature. The 
calculation is based on the variational principle in the form of a perturbative series expansion for 
p(T). A peculiar behavior of p(T) in L1V2O4 is related to properties of low-energy spin fluctuations 
whose T-dependence is obtained from SCR theory. 

PACS numbers: 71.27.+a, 71.10.-w, 72.10.Di, 72.15.-v 

I. INTRODUCTION 

The metallic vanadium oxide Li V2O4 has attracted much attention after a heavy fermion behavior in this 3d-electron 
system was discovered^— The cubic spinel LiV204 has the pyrochlore lattice of vanadium ions (in the mixed valence 
state V 3,5+ ) and shows metallic conduction and no long-range magnetic ordering for any measured temperatures 
at ambient pressure. So far the origin of the heavy fermion quasiparticle formation observed in this compound 
for T < 10K remains to be a controversial subject, however, effects of electronic correlations and the geometrical 
frustration of the pyrochlore lattice are supposed to be key aspects of the problem. 

The quasiparticle mass enhancement is expected when a metallic system is driven by strong electron correlations to 
a vicinity of a charge and/or spin phase transition at low T. In that case, the charge/spin disordered ground state on 
the metallic side of the transition in the strongly correlated system LiV204 is sustained because a long-range order 
with a particular ordering (critical) wave vector Q c is prevented by the geometrical frustration. Expressed differently, 
the system cannot choose a unique wave vector of an ordered structure which minimizes the free energy. Instead, it 
is frustrated between different structures with different critical wave vectors Q c 's and equally low free energy. For 
instance, low energy spin fluctuations are expected to be present in a very large region of momentum space which is 
the signature of frustrated itinerant magnetism. This is in contrast to non-frustrated systems where the fluctuations 
are confined to the immediate vicinity of a unique incipient ordering vector. 

This scenario for frustrated itinerant magnetism was recently investigated in detail for LiV204 by present authors. 4 . 
An analysis of inelastic neutron scattering (INS) measurements 5-7 and calculations of the dynamic spin susceptibility 
allowed us to suggest the location of the paramagnetic spinel L1V2O4 close to a magnetic instability. This was 
achieved by developing the RPA theory of spin fluctuations based on ab-initio band structure calculations and an 
on-site Coulomb interaction of 3d electrons. Close to the critical value of the interaction strength, low energy spin 
fluctuations develop throughout a large shell in momentum space. They may be mapped to an effective low energy 
paramagnon model which describe low-temperature INS results 5 - - — accurately. From the comparison at T — >• 0, the 
parameters of the model (peak energy, weight and extension in momentum space) are fixed. Using the selfconsistent 
renormalization (SCR) theory^ - — which includes mode-coupling of spin fluctuations, the finite T properties of INS 
spectral shapes, uniform and staggered susceptibility, as well as NMR relaxation rate, have been explained 4 -!^. From 
this analysis we concluded that L1V2O4 can be regarded as a nearly antiferromagnetic (AFM) metal and its unusual 
low-T properties have to be related to a peculiar structure of the paramagnetic ground state with strongly degenerate 
low-energy (slow) AFM spin fluctuations. 

In the present study, our main concern is to explain the low-temperature, T < 40K, electrical resistivity p(T) 
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measured on single crystals of L1V2O4 and reported by Takagi et alM* and Urano et alJ^. A Fermi-liquid behavior 
p(T) ~ T 2 for T < 2K and a more slow increase of p{T) for higher temperatures were found. Measurement o 10 ' 11 
revealed a noticeable change in physical properties of LiV204 for T > 40K, including a Curie- Weiss magnetic suscepti- 
bility x{T) and a highly incoherent transport, which is, however, beyond the scope of present theory. In our approach 
we will use the effective low energy paramagnon model for spin fluctuations whose parameters are completely fixed by 
the comparison with INS. Only two more pheonomenological parameters characterising the impurity and paramagnon 
scattering mechanism will be needed. 

From an analysis of the optical reflectivity and conductivity measurements, Jonsson et alr^ and Irizawa et al^ 
inferred that the conducting electron system in LiV204 at ambient pressure is located close to a correlation-driven 
insulating state. Under the applied external pressure^ the system undergoes a metal-insulator transition accompanied 
with a charge ordering and a structural lattice distortion in the insulating phase. The observed complicated phase 
transformation display properties different from those expected from a first-order phase transition or a conventional 
metal-insulator transition. We note that because of quarter- filling of the electronic t^ g bands in L1V2O4, electron 
correlations due to inter-site Coulomb repulsion have to play an essential role in the observed transition.— Under such 
conditions, a microscopic mechanism for the heavy quasiparticle formation on the metallic side of the transition has 
to be clarified. However, in the present analysis which is concerned with resistivity under ambient pressure the effect 
of inter-site Coulomb interaction and a slowing-down of charge fluctuations will not be included. 

It was realized some time ago that the calculation of low-T transport properties in nearly AFM metals is a rather sub- 
tle issued - — In clean systems with a peculiar AFM ordering vector Qafm the quasiparticle scattering by quantum- 
critical spin fluctuations is strongly anisotropic. The strongest scattering occurs near the "hot spots" of the Fermi 
surface (FS) connected by Qafm and the main contribution to the electrical conductivity is due to quasiparticles 
from the "cold region" of the FS. In that case, if the system is at some distance from the AFM quantum critical 
point, the low-T scattering rates are proportional to T 2 and the Fermi-behavior p(T) — pi mp ~ T 2 is realized. As was 
first pointed out by Rosch ) 16 i if a small amount of disorder is present, an interplay of strongly anisotropic scattering 
due to critical spin fluctuations and an isotropic impurity scattering may complicate the picture producing several 
different regimes for p(T) — pi mp ~ T a with the exponent between 1 < a < 2 in the low-T region. 

To describe temperature-dependent electrical resistivity in L1V2O4, we suggest that quasiparticles are scattered by 
AFM spin fluctuations almost equally strongly over the FS and effects of anisotropy are weak. This is related to a 
peculiar distribution of dominant AFM spin fluctuations with ordering vectors Q c 's forming a largely isotropic dense 
manifold in k space in this compound. Low-temperature evolution of interacting spin fluctuations in L1V2O4 can be 
successfully described within the SCR formalism as presented in our previous studies^ As explained in Refill, the 
set of model parameters of the SCR theory were obtained^ from neutron scattering data£~— and used to describe^ the 
temperature and pressure evolution of the spin- relaxation rate 1/TjT observed in the NMR measurement^! of the 
low-T spin dynamics in L1V2O4. In this work, the theory is extended and applied to give an explanation of the low-T 
electrical resistivity p{T) in this compoundi 10 ' 11 



As follows from experimental observations^"— ' 10 ' 1 1 the concept of the Fermi quasiparticles for charge carriers in 
the metallic spinel LiV2 04 is valid for sufficiently low temperatures, T < 30 K. In this regime, we assume that the 
dominant scattering processes are given by low-energy AFM spin fluctuations and impurities. In the linear response 
theory, in an applied electric field E the quasiparticle distribution function / k is linearized around the equilibrium 
Fermi distribution according to / k = — ^^df^/de^. The electronic transport can be found from the Boltzmann 
equation 



The scattering operator Wkk' can be expressed through the total equilibrium transition probability "Pkk' = V^f 

as (&b = h = 1): 



II. VARIATIONAL PRINCIPLE FOR p(T): GENERAL CONSIDERATION. 




(1) 




(2) 



provided the spin fluctuations are in thermal equilibrium, i.e., there is no drag effect. 
For the elastic impurity scattering one has 



P£* = 2™ 4 |T kk ,| 2 <5(e k _ £k ,)/0(i _ /o ). 



(3) 
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To a sufficiently good approximation, the T-matrix in Eq.® is frequently assumed to be a constant Xkk'| 2 ~ Vi mp 
and fiiVi mp , where rij is the impurity density, is regarded as a free parameter to be chosen so as to give a realistic 

value of the measured residual resistivity pi mv - We avoid this approximation and treat below matrix elements of V^f 
generally. 

For the spin-fluctuation (s/) scattering one ha o 15 i 17 ' 22 ' 23 

Kl> = 3^ 2 //k(l - /kOK^k - ekO + l]Im X (k - k', e k - e v ) , (4) 

where n(e) is the Bose distribution function, x (q, e) is the dynamical spin susceptibility describing the low-T para- 
magnetic state of L1V2O4 and J s f is an effective coupling constant which is the second free parameter. It is worth 
emphasizing that in the present study the other parameters of the phenomenological SCR theory determining the 
behavior of x (Qi e ) are considered to be known and fixed from a fit to the data of inelastic neutron scattering 
measurement— on L1V2O4, as discussed in RefH. 

Following the standard notation 2 ^, the Boltzmann equation ([1]) can be rewritten in the form Xk = Wkk'^k'- 
Then, the electrical resistivity can be obtained by minimizing a functional 2 ^ 



\(^X(E = 1W 



(5) 



Here, E = 1 means the unit electrical field and the scalar product of two functions $k and ^k is defined as ($, = 
Sk ^k^k- In fact, in Eq. (|5j) the k-integration over the actual FS is implied, which follows from the property of the 
scattering operator Wkk' and the explicit form of A"k — e(Evk)(— df^/de^). 

A way to search for a variational solution of Eq.Q for the deviation function $k is to expand it in a set of the 
Fermi-surface harmonics (FSH) 0L(k): 

$k = 5>L<Mk), (6) 

L 

where tjl are variational parameters and L is a convenient composite label that includes numbering of different sheets 
of the FS in LfV^CU. The FSH's are defined 25 ' 26 as polynomials of the Fermi- velocity Cartesian components v£. That 
is, for each integer N > one has to construct (N + l)(N + 2)/2 polynomials (v£) 1 (v^) m (v^.) n with l,m,n > and 
I + m + n = N, and orthonormalize them on the actual FS, (4>l', 0l) — ^L'L- The resulting polynomials forming a 
complete set {</>i(k)} of basis functions are classified according to different irreducible representations T of the lattice 
symmetry point group. In general, for a given T there are subsets of different functions, {<j) r L (k)} , {<j) s L (k)}, etc., which 
transform according to the same T. Then, for any pair of partner functions, 0£,(k) and 0^„(k), belonging to different 
subsets, but transforming according to the same row of T, one has ((jf L ,,W(j) s L „) ^ 0. The other off-diagonal matrix 
elements of the scattering operator PFkk', including those connecting different irreducible representations, vanish by 
symmetry arguments and, hence, the scattering operator has a block-diagonal form (see the discussion by Allen— and 
references therein). 

A minimum of is achieved in the class of odd functions, <£t(— k) = — 0i(k); only these basis functions are 
included in the expansion (J6j) . Recalling the cubic symmetry of the LiV2 04 lattice structure, we assume without 
loss of generality that the applied electric field points in the x direction, which immediately distinguishes one of the 
first-order FSH's: <j) lx (k) = v£/ («) 2 } 1/2 , where («) 2 ) 1/2 is for the root -mean-square on the Fermi surface. 

A general strategy in describing the physical resistivity p(T) as a solution of the variational equation ([5]) in most 
of the metallic systems, including those with complicated electronic band structure, is to truncate the expansion ([6]) 
by keeping in it only a few of FSH's. Following the common practice, one may start the analysis with the lowest, 
first-order variational solution, 4>| c °' ) ~ cf>i x (k) , which is a fairly good approximation provided the anisotropic effects of 

the quasiparticle scattering are weak. As usual, here the anisotropy of the scattering operator Wkk' = W^£? + W^, 
means that the transition probabilities depend not only on the mutual angle between the momenta k and k', but also 
on their position with respect to the crystallographic axes. Anisotropic effects, as well as a complexity of the actual 
FS, can be partially caught in the calculations by keeping in the expansion (O a selected number of higher-order 
FSH's. 

To go beyond the lowest-order solution for in the simplest manner, the following approximate assumption 
can be made: the off-diagonal matrix elements Wll' — {4>L,W(f>L<) are small compared to diagonal ones, Wll and 
Wl'L'- Then, the variational solution to Eq.©, being written in the familiar form 24 as p^ 1 = Xf(<fii x , W -1 </>i x ), 
where X\ = (4>i x , X), can be expanded in a perturbation series^ 6 - 



x\ 



st*' W 1x ,lWl,i x st-'W^lWluWl' 

2—1 W, Wr r 



L W lx , lx W LL f£ W 1x , 1x WllW l >l> 
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where the primes on the sums means that the terms with L, V = lx and L = L' are excluded. 

For T ->• 0, from Eqs.© and (0} one has W s L f L , -> 0, both for L = V and L ^ L', and the Eq.© reduces to 

p(T 0) = p 4mp « XfV^LC (8) 

where the constant £ stands for brackets in Eq.0 with VFll' = W™l P . Its value is less than unity, < £ < 1, since 
the inclusion of higher-order terms leads a lower estimate of the upper bound for p. The expression given by Eq.® 
approximates the experimental valu o 10 ^ 11 of the residual resistivity ~ 32 pilcm in a low-T fit procedure. 

For the further purposes, we note that relations between W™\ x and the other diagonal matrix elements W^™ p 
cannot be generally established. In particular, a strong inequality W^2 P ^> W^\ x for some L ^ lx is not excluded, 
which does not invalidate our previous assertions. Actually, irrespective of a relation between W™\ x and W^™ p , the 
series expansion ([7} starts with the matrix element W™ p x due to the requirement that the variational solution for 
p is given by the diagonal matrix element of the inverse scattering operator W~ x between the same first-order FSH, 
i.e., p- 1 W-^ix). 



III. ELECTRON SCATTERING BY SPIN FLUCTUATIONS IN LIV 2 4 . 

An open question is: Whether one may rely an analysis of the physical resistivity on the series expansion ([7J for 
T > 0? Since the impurity scattering is thought to be highly isotropic, one of the underlying assumptions that 
W^2? /W^2 P <C 1 for L' ^ L, has to be fulfilled. Here, the appearance of some off-diagonal matrix elements W^2? can 
be explained mostly due to a complex character of the multi-sheet FS in L1V2O4. Below we examine how properties 
of AFM spin fluctuations are related to those of the spin- fluctuation scattering operator W^, in LiV204, and show 

that the smallness of the off-diagonal elements W^J L , with respect to diagonal ones, W S L { and W S J V , seems to be a 
plausible assumption as well. 

From Eqs.Q and ((4]), any diagonal or off-diagonal matrix element W^ L , allowed by symmetry arguments can be 
written as follows 

w ll> = -jr E t^( k ) - ^( k ')l v & [<Mk) - <Mk')] . (9) 

kk' 

The use of the definition (|4]) for V^, leads to 

w ™ = h ( h) 6 / de f $ j ^7 {<i>L (k) - ^ L (k ' )] psf (k ' - k; e) [0L ' (k) - 0L ' (k ' )] ' (10) 

where the standard replacement — > f de <j> ePfc/[(27r) 3 Uk] for a unit volume together with the relation /°(ej c )[l — 

/°(ek')][n(ek - £k') + 1] = [/°(ekO - /°(ck)]n(ek' - ek)[n(e k ' - e k ) + 1], and the approximation [/°(e) - /°(e')] « 
(e' — e)(—df°/de) are used; two-dimensional integrations over fc and k' are restricted to the FS. Then the kernel 
V(k' - k, e) in Ea.([T0]l takes the form 

V sf {k' -k,e) = 3J s 2 /e n(e)[n(e) + l]Imx(k'-k,e). (11) 

This is the well known form of conduction electron scattering from spin fluctuations. The latter are enhanced by the 
nearly critical Coulomb interaction of 3d electrons, which leads^. to a paramagnon expression for \ s -' ~ k, e) whose 
parameters are fixed from INS results at T — > 0. At low temperatures, the low-energy (e ~ 1 meV) dynamic spin 
susceptibility x(q'i £ ) i n L1V2O4 shows maxima around the critical wave vectors q' = Q c forming a rather dense 
manifold {Q c } in k space.— To take into account explicitly all scattering processes due to dominant spin fluctuations, 
it is helpful to make in Eq. trTTI) the following substitution 

Im X (k' - k, e) ~ T I 7?T3 5 ( k ' k Q<= q)Im X (Qc + q, e) , (12) 

which ascribes particular weights Im% (Q c + q, e) to the quasiparticle scattering processes whose wave vectors, k' and 
k, at the Fermi surface satisfy the relation k' — k = Q c + q. In Eq. (fT2|) . the summation is over the entire set {Q c } 
of the critical wave vectors and their neighborhoods, |q| <C |Q C |- In total, this involves a broad region in k space, 



5 



a 



0.15 



0.10 



0.05 



0.00 




y Q (0) - 0.044 
bx c = 0.01 
T = 60 (K) 
T A = 220 (K) 



0.5 
t 



1.0 



FIG. 1: The solution of the basic equation of the phenomenological SCR theory for the reduced inverse static spin susceptibility 
UQ(t) as a function of the reduced temperature t = T /To; phenomenological parameters required to obtain j/g(i) in LiV2 04 are 
given in the text. 



where the dominant AFM spin fluctuations are distributed, and Imx (Q c + q 7 e) in Eq. (fT2l) does not much depend on 
a direction of Q c . The resulting distribution differs strongly from that occurring at low T in most of nearly AFM 
metals where the low energy susceptibility \ (q, e) is usually peaked around a discrete ordering wave vector Qafm- 
The use of the above arguments allows us to write down the matrix element W^f L , in a factorized form 

W S L {, *C W F{T), (13) 

where 

Cw = (h) f ? / [Mk) ~ Mk ' )] M ™ [<Mk) ~ 0i ' (k ' )] ' (14) 



HT) = \\^\ (0 m ( £ ) H £ ) + ^ lm X (Qc + q, e) , (15) 
o 

and the matrix M^, is denned as M^, ~ 3J^- X){q } k Qc)- The matrix is invariant under simultaneous 

operations of the lattice point group on both k and k', since the manifold {Q c } is an invariant as well. The wave 
vectors Q c are along all high symmetry directions in the Brillouin zone (BZ) and their end points are lying on a 
closed surface of a mean radius |Q C | ~ 0.6A -1 which is referred to as the "critical" surface. 4 The factorization of 
matrix elements introduced by Eq. (|13|) implies that the quasiparticle scattering by spin fluctuations with different Q c 
at the "critical" surface provide nearly identical contributions and, therefore, only one representative wave vector Q c 
appears in Eq. ([15[) . 

As discussed in RefS the high directional degeneracy of Q c 's results from complexity of the electronic band 
structure and the geometrical frustration of the pyrochlore lattice structure of LiV 2 04. In this respect, the low- 
T spin-fluctuation scattering mechanismi^r— operating with a peculiar ordering wave vector Qafm differs from 
that occurring in the paramagnetic spinel L1V2O4. In the former case, the quasiparticle scattering is a strongly 
anisotropic one, leading, for instance, to "hot spots" at the Fermi surface. Instead, the above analysis suggests that 
the quasiparticle scattering by spin fluctuations in L1V2O4 is largely isotropic one. In that case, the diagonal matrix 
element W%f L prevail over the off-diagonal ones. This fact and similar arguments mentioned above for the impurity 
scattering justify the applicability of the perturbation series expansion, Eq.(J7]), for p(T). 

IV. CALCULATION OF p(T) BASED ON THE SCR THEORY OF SPIN FLUCTUATIONS IN LIV2O4. 

First, based on the expansion Q, the resistivity p(T) can be expressed as follows 

p(T)- Pimp = pW(T) + Ap(T), (16) 
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where the spin-fluctuation contribution • (T) = X 1 2 W^ lx represents the lowest-order solution and Ap(T) is a 
correction due to higher-order terms in Eq. (|7|) . Their variations with T can be found by calculating the function J-(T) , 
Eq. (|15[) . entering the matrix elements of the spin- fluctuation scattering operator, Eq. (|13[) . Although the subsequent 
derivation of an explicit form of J-(T) has much in common with earlier studie d 23 ' 27 of nearly AFM metals, essential 
features specific to L1V2O4 have to be emphasized. In particular, within the SCR theory of spin fluctuations the 
imaginary part of the dynamic spin susceptibility can be parametrized as follows^ 



Im X (Q c + q,e;T) = 



1 



4ttT a T 



VQ (T) + (qW/q B ) 2 + b(^/q B y 



(e/27rToY 



(17) 



Here, for a given Q c , qW and are the components of q parallel and perpendicular to Q c , respectively; qs is the 

1/3 

effective radius of the BZ boundary given in terms of a primitive cell volume vo as qs = (67r 2 /«o) • The parameters 
Ta — 220K and To ~ 60K characterize the widths of the momentum and energy distributions of spin fluctuations, 
respectively; a small parameter b takes care about a strong anisotropy of the distribution in k space. Next, the 
reduced inverse susceptibility at Q c is defined as j/q (T) = \2Tax{Qc',T)\~ 1 assuming a "spherical" approximation, 
i.e. x (Qc! T) does not depend on a direction of Q c . 
With the insertion of Eq. ([T7|) into Eq.(TT5j). we get 



HT) = 



27rTo f d 3 q 



d\n(\)[n(\) + 1] 



2^T n 



A 2 + (27ni q ) 2 T A J (2tt) 



d 3 c 



(18) 



Here, the last equality defines I(uq), where 



Uq(t) 



VQ (t) + (q ll /q B ) +b(q L /q B ) 



(19) 



and t = T/Tq is the reduced temperature. At the next step, the function I(iiq) can be expressed as 



no 



i q 2"U 2 



(20) 



where ij){u) 1S t ne digamma function and ip'( u ) = dtp(u)/du. 

The integration over q in Eq.(fT8"l) are performed by using the same prescriptions as in Ref|8|, which yields 
(2ir)~ 3 J d 3 ql(uq) — c/(T/Tq), where c is the known dimensionless factor, c » 4, and^ 7 - 



f(T/T ) = — I dz 



bx c 1 , 

In 

t 2 



y Q {t) +z 2 + bx c 
VQ(t) + z 2 



Vq{t) + z 2 + bx c ( y Q (t) + z 2 + bx c \ _ y Q (t) + z 2 ( y Q (t) + z 2 



lnT 



y Q (t) + z 2 + bx c 



t ) t 

VQ{t) + z 2 ^ 



Mr 



t 



(21) 



where T(u) is the gamma function, the cutoff z c ~l/2 and the remaining parameter bx c ~ 10 2 . 

For a given UQ{t), the expression (|2ip. as a function of temperature, can be calculated numerically. This determines, 
according to Eqs. (fr^ . (rT8|) - (l2T1) . an evolution with T of any non- vanishing matrix element W S J L , in the whole range 
T < 40K, where the SCR theory for the AFM spin fluctuations in L1V2O4 is proved to be valid&. Here we utilize 
the known functional form for yQ(t) 7 Fig.l, obtained by solving the basic equation of the SCR theory developed in 
RefH to explain results of inelastic neutron scattering measurements on L1V2O4. The solution shows that 2/q(t) is a 
monotonically increasing function of temperature; the limiting value j/q (i — > 0), was found to be Vq(0) =0.044. The 
use of the energy scale, T* = 2irT yq(Q) w 16K, which is the relaxation rate of the low-energy spin fluctuations, 8 
is helpful in recognizing two regimes with different power-law behavior of j/q (t) . Actually, for T <C T* , one obtains 
[2/q(*) ~ 2/q(0)]/?/q(0) < 1, which leads to the quadratic behavior of f(T/T ) — T 2 . For T* < T, a smooth, nearly 
linear, t-dependence of yqit) results in a peculiar monotonic temperature increase of /(T/Tq), as indicated below. 
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FIG. 2: T-dependence of /(T/To) defined by Eg. (|2f [I and calculated with the use of the reduced inverse susceptibility uq(T/To) 
shown in Fig.l. 



First, we analyze the lowest-order approximation to the low temperature (T < T*) resistivity, 



P(T) W Pim V + Psf(T), 



where 



$(T) = A$f(T/T ) 



(22) 
(23) 



and compare its T-dependence with that of the observed^ experimental resistivity. Hereby = (2ttcT /Ta)X 1 2 C\ Xt \ x 
is an adjustable parameter (together with pi mp ) in a low-T fit procedure using the calculated /(T/T ) shown in Fig. 2. 

One may see that the function /(T/To) nearly precisely follows the quadratic dependence, /(T/T ) = c\T 2 with 
ci = 0.0033, for T < 2K < T*, where the Fermi-liquid behavior [p exp (T) — p imp ] = AT 2 in LiV 2 4 was reported ,±&±i 

From the low-T fit procedure, as shown in Fig. 3, the parameter A^) was found to be A^) = 666.7 pQcm, which 
corresponds to the observed coefficient of the T 2 term, A = 2.2 pflcm/K 2 . 

With increasing T and starting from T i=a 2K, both the calculated p^j '(T) and the measured resistivity [p exp (T) — 
Pimp] show gradual deviations from the T 2 behavior, however, with somewhat different rates. Specifically, starting 
from T » 2K one obtains the growing discrepancy [p exp (T) — pi mp ] — p^)(T) = Ap{T) < 0, which indicate that the 
higher-order corrections involved in Ap(T) cannot be longer neglected. Remarkably, a negative correction, Ap(T) < 0, 
is consistent with the variational principle requiring that an extension of the involved basis functions should lead to an 

(2) 

improved upper bound on p(T). Actually, let us consider the lowest-order correction Ap s J(T) which can be written 
from Eq.flUl as 



A^(T) = -^' 



X 2 ^ 



x,L + vv lx,L) 



w. 



sf 



LL 



wTT 



(24) 



where the second term in brackets, being already involved in Pi mp , is now subtracted to ensure that Apy/ (T -> 0) = 

0. Note that a T-dependence in the right-hand side of Eq. (|24"|) is entirely due to W^f L , — Cll' J~(T), both for 
L' = L and L' ^ L. For sufficiently low T, a denominator in the right-hand side of Eq. (|2"4")l can be approxi- 
mated assuming that CllJ~(T) <C W^2 p and T(T) ~ T 2 , which immediately leads to small leading correction, 
a p1/( t ) ~ aT2 - \b\ Ti + 0(T 6 ). Here, the first quadratic term can be adopted by changing slightly a value of the 
adjustable parameter A *), while the next term, — |£>|T 4 , provides the required negative correction to the first-order 
result p$(T). 



An extension of the above analysis to higher temperature could be possible if one establishes reliable relations 
between numerous matrix elements involved in Eq.([7|). The following plausible assumption can be made based on the 
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FIG. 3: Theoretical fits to the experimental data (Reffiol. open circles) for the electrical resistivity in LiV204 in different 
temperature regions. In the low-T limit, the resistivity is described by Eqs, (l22[) . (|23l) with pi mp = 32(/if2cm) and = 
666.7(/iSlcm). In the high temperature region, 20 K < T < 40 K, where the SCR theory is still valid, the resistivity is 
approximated by Eq. (|25[) with pi mp + B = I17(/if2cm) and A s f — 280(/^f2cm). 



fact that T{T) = (2ttcT /Ta)/(T '/T ) is a rapidly growing function of T, Fig.2. For instance, F(T ~ 30K)/J"(T ~ 
IK) ~ 10 2 . We suggest that the limit, CllJ~(T) ^> W L T , can be achieved at T < 40K, i.e near the border where the 
SCR theory of spin fluctuations in L1V2O4 is still valid. With this assumption, one obtains, for instance, the following 
estimate for the second order correction, Apf»(T) w — [^^^^{Cix^Y/CLL^iT) + const. Then, by doing in 
the same manner and after collecting all leading terms in the expansion ([7]), the physical resistivity in L1V2O4 for 
2K<C T < 40K can be approximated by the following simple functional form 



p(T)^p imp + B + A s ff(T/T ), 



where 
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LL 
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LL 1 



'Ci x lCll>C_ 



L',lx 



CrrC 



LL^L'L' 



Cix.L 

Cll 



1 



Clx, 



Cll 



(25) 



(26) 



(27) 



By noting a close similarity between Eq. (|23|) and Eq. (|2"21 , it is worth emphasizing completely different T-dependence 
of f(T/To) in the low- and a high-temperature regimes. Moreover, the factor A s f is subjected to a special constraint 



with respect to A *). Actually, for T < 2K <C T* the first-order term in the series (|26|) is only needed, A 



1(1) 
<A 



(2ttcT /T a )X^ 2 Cix,ix, while for 2K < T < 40K, the factor A s f is given by the full series (|2"6"j) and, hence, A s f ^ ^ s f 
is required. An estimate for a shift B in Eq. ([25|) . which is the second adjustable parameter in a high-T fit procedure, is 
discussed bellow. In Fig. 3, the physical resistivity [p exp {T] — pi mp ] is compared with the predicted behavior, Eq. ([2~5|) . 
for T > T*(« 16K). A satisfactory coincidence between the experimental data and calculated results is achieved 
for 20K < T < 40K with two fit parameters A s f = 280u£lcm and B = 85uilcm. While the expected constraint, 
A s f < A^f is fulfilled, the obtained large value of B ~ 3pi mp means that the Matthiessen rule 2 ^ is severely violated. 
A possible mechanism for this effect is the following. Recalling the estimate, Eq. (l27j) for B, together with relations 
and [Cix.l/Cll) -C 1, we suggest that for some L ^ Ix one has W L T ^ WiTix' 



Pimp ~ Wi x ix and (Cix,l/Cll) -C 1, we suggest that for some L ^ lx one has W LL ^ ^lx.ix' which explains 
why an estimate B ~ pi mv is feasible. 

So far the special attention has been paid to two limiting regimes of low- and comparatively high-temperatures 
(T < 40K), where the series expansion, Eq.©, for p(T) reduces to very similar forms, Eqs. (|2"2"|) and (|23|) . requiring 
two adjustable parameters for a fit procedure in each regimes . We insist that Eq.© should provide the interpolation 
T-dependent function between the low- and high-T limits as well. However, for intermediate temperatures, one has 
CllF{T) ~ M / 2™ p , and the corresponding fit procedure, though being possible, would require a larger number of 
adjustable parameters. This could hardly give more insight into the problem under consideration and therefore we 
omit such a procedure in our discussion. 
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V. CONCLUSION 



We have calculated the electrical resistivity p(T) in the paramagnetic metallic spinel L1V2O4 treated as a nearly 
AFM Fermi-liquid for temperatures T < 40K . Impurities and strongly degenerate temperature-induced low-energy 
AFM spin fluctuations were supposed to provide two main sources of the quasiparticle scattering and the resistivity. 
The self-consistent renormalization theory developed earlier was applied to derive explicitly the temperature-dependent 
matrix elements of the spin fluctuation scattering operator. The absence of hot spots of the Fermi surface and a largely 
isotropic character of the quasiparticle scattering was deduced from a peculiar, nearly spherical, shape of the spin- 
fluctuation distribution in the momentum space for the paramagnetic ground state in LiV2 04. Comparatively weak 
anisotropic effects were assumed to originate mainly from a complex many-sheet structure of the Fermi surface in 
this compound. The assumption allowed us to use the variational solution for the Boltzmann equation in the form 
of a perturbative series expansion for p(T). Our theory remains to be a phenomenological one since unknown model 
parameters were found from the best overall fit to the temperature-dependent p exp (T) measured on a single crystal 
ofLiV 2 4 . 

The resulting theoretical expression for p(T) was shown to take very similar simple forms in two limiting regimes 
for spin fluctuations, which describe successfully experimental results for p(T) with a minimal set of two adjustable 
parameters in each regime. These include the low temperatures, T <C T* (where T* « 16K is the characteristic energy 
scale of spin fluctuations), and somewhat higher temperatures, T* < T < 40K, respectively. 

For T > 40K the SCR theory of AFM spin fluctuations in L1V2O4 is no longer valid. As discussed in Ref.S, and 
evidenced from experiment)* 7 ' 10 ' , with increasing T the AFM fluctuations at |q| ~ Q c are suppressed and no more 
distinguished from those at other wave vectors in the BZ; the system enters a spin localized regime compatible with 
the Curie-Weiss behavior of x(q = 0) observed in LiV2C>4 for T > 50 K. An explanation of incoherent transport 
properties in this regime remains to be a challenging problem. 
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